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FORTRAN IV PROGRAMS 
FOR SUMMARIZATION AND ANALYSIS OF 
FRACTURE TRACE AND LINEAMENT PATTERNS 


ABSTRACT 


Systematic and detailed analysis of lineament and fracture 
trace patterns has long been neglected because of the 
large number of observations involved in such an analysis. 
Three FORTRAN IV programs were written to facilitate 
this manipulation. TRANSFORM converts the initial 
fracture map data into a format compatible with AZMAP, 
whose options allow repetitive manipulation of the data 
for optimization of the analysis. ROSE creates rose 
diagrams of the fracture patterns suitable for map 
overlays and tectonic interpretation. Examples are 
given and further analysis techniques using output from 
these programs are discussed. 
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FORTRAN IV PROGRAMS 
FOR SUMMARIZATION AND ANALYSIS OF 
FRACTURE TRACE AND LINEAMENT PATTERNS 

INTRODUCTION 

Photogeologic fracture traces and lineaments^ or other closely related syno- 
nyms (micro- and macrofractures^, megajoints^, lineaments'^, etc. ) have 
enjoyed a resurgence of interest in the last several years. They have been 
noted on both large and small scaled imagery of the Earth and surrounding 
planets^, and in places on Earth, can in many cases be seen throu^ consider- 
able overburden®. Their occurrences and associations with the local and 
regional tectonic stress patterns are often enigmatic. 

One of the major problems of fracture analysis is in handling the large number 
of observations. Not only does this problem render repetitive manipulation 
impractical, but also makes the varying of classification parameters to 
achieve an efficient and sufficient classification cumbersome. As an example, 
if one wishes to analyze the data for local perturbations in a fracture pattern, 
the smallest grid size (siunmation area) that would achieve statistical 
stability might be desired. Relatively small areas (200 square miles) may 
contain several thousand fractures. This paper will present a series of 
computer algorithms which will greatly facilitate fracture (fracture trace) 
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analysis and will present some results of one particular study. It should 
be emphasized that the same techniques may be applied to lineament analysis^ 
bearing in mind the changes in scale and their possible significance. 

FORTRA N IV PROGRAMS 

Purpose 

Three programs were written in FORTRAN IV. Their major purpose is to 
process and summarize in various fashions the large number of primary 
observations of fracture trace orientation, length and position, while varying 
the classification parameters so that further statistical or structural analysis 
techniques might be applied. Input to this group of programs consists of 
fracture pattern maps plotted on a suitable base, while output consists of 
frequency-azimuth histograms and/or rose diagrams for the mapped area. 
Small areas (subsets) of the original mapped areas may be examined by a 
gridding technique and additional statistical parameters generated. 

TRANSFORM Program 

TRANSFORM performs the initial data treatment, converting the fracture 
traces on the base map into a format acceptable for the AZMAP program. 

A Cartesian coordinate system is established with its origin in the upper left- 
hand comer of the map. The X-axis is latitudinal and positive to the right 
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while the Y-axis is meridional and positive downward. North is assumed at 
the top of the map. Due to the convergence of meridional lines, there will be 
some discrepancy between true north and the Y-axis for any given point on the 
map, however, this variation is limited to a total of 2 degrees on a 1:250000 
Merpator plot map at 45 degrees latitude. 

With the establishment of this coordinate system, the begiiming and end of each 
fracture may now be indexed. Digitizing equipment, which automatically 
records coordinates either on standard Hollerith computer cards or magnetic 
tape, is used to reduce the fracture map to a form usable by the program. 

The digitizer "bull’s eye" is set at the origin and a scale is established 
accordii^ to the requirements of the individual machine. Once this is done, 
the beginning and end points of each fracture may be referenced by placing 
the "bull’s eye" on the beginning (greater Y value) point, recording the value 
as per the digitizer format on Hollerith cards or magnetic tape, followed by 
the end point (smaller Y value), repeating the recording procedure. In order 
to allow the operator to digitize the fractures in a logical manner, it is 
suggested that each fracture be given an inventory number prior to digitiza- 
tion, beginning with one and consecutively numbered. Most digitizers will 
generate consecutive numbers on a display panel, thus allowing the operator 
to constantly, check if any fractures have been missed in the digitizing 
procedure. 
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The program treats each fracture as a vector in map space and generates 
parameters used in classification techniques of the AZMAP program. This 
procedure is only required once for each digitized map, and subsequent 
repetitive treatments are done in AZMAP and ROSE. The Appendix contains a 
program listing, abstract, and detailed instructions for data processing. Less 
than 30 seconds were required to process 1500 vectors on an IBM 360/67 
computer. 

AZMAP Program 

AZMAP uses data generated by the TRANSFORM program. In a mechanical 
sense, the program places an orthogonal grid over the map parallel to the 
designated X- and Y- map axes. The operator selects the summation area size 
(cell) over which the fractures will be summarized by specifying the X- and Y- 
axis grid cell size. The computer then scans all the fractures, determining 
whether each fracture falls within the grid cell, incrementing the cell by 
operator- supplied values in both the X and Y directions, until the total 
designated area has been covered. By incrementing at submultiples of the 
cell size a sliding average technique may be employed. 

Either of two summary techniques may be used. Subroutine MID counts the 
whole fracture as falling within the cell if its midpoint falls within that cell, 
while Subroutine PART considers only that portion of its length which lies 
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within the cell. The choice of either subroutine depends iqjon grid cell size, 
size of the linear features mapped, and goals of the operator or experiment. 

The operator also selects the number of azimuth classes into which the data 
will be summarized. Up to 90 classes may be specified, their values being 
incremented from 270 through 0 to 90 degrees. Thus, if a fracture trace lies 
within a specific cell, it is then added to the appropriate azimuth class within 
that cell according to the subroutine selected, and a histogram is plotted. 

Data is summarized both as density (total length of fractures within each cell) 
and frequency (number of fractures within a cell). Depending upon the tech- 
nique used, either the whole length (MID) or the portion (PART) is summed 
for density while a whole unit is added to frequency in either case. Therefore, 
if a fracture trace extends into or through three cells, frequency will be 
incremented by one in all three cells if subroutine PART is used, 

A Chi Square (y^) test is performed on the resultant frequency-azimuth histo- 
gram for each grid cell if so desired, testing the distribution for randomness 
(a rectangular distribution, where all classes have an equal chance of 
occurring)^. Care should be exercised in the interpretation of the y2 test, 
because the lower limit of reliability is reached when the expected frequency 
for each azimuth class is less than 2 units where: 

Total number of units 

Expected Frequency = — 

Number of azimuth classes 
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Caution is also advised for the density test because the test results will be 
dependent on the scale of the units chosen in this classification. 


Figure 1 contains an example of the line printer output for one grid cell. In 
addition, punched card or magnetic tape output may be generated for use 
either in the ROSE program or additional statistical programs that will be 
described but not documented in this paper. The Appendix contains a source 
listing, abstract, and detailed instructions on the program's use. Timing 
considerations are dependent upon the svimmary technique. MID is more 
efficient and requires less than 60 seconds for summary of 1500 fracture 
traces into 49 grid cells, while PART requires about 90 seconds on an 
IBM 360/67 computer. 

Although AZMAP requires a relatively large computer, the program may 
be modified to fit on smaller computers by reducii:^ the size of the arrays 
designated for storage of the fracture data (dimensioned at 2000) without a 
resultant decrease in program efficiency. If time is of no objection, then 
the above mentioned arrays may be removed and with a little rewriting, the 
fracture data could be stored on tape and repetitively scanned and rewound 
for each grid ceR summarized. 
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TEXAS STUDY AREA* 1/16TH AREAS* 1/2 CELL INCREMENT 

EACH G^IO C£LL is 240 MM. ( 3.576 MILE< 5 )> BY 260 MM. f 3 .B 74 M 1 LE«S)| 

PROGRAM USES SU3R3UTINE MID* CONSIDERS hfHOLE VECTOR AS BEING WITHIN CELL IP ITS MIDPOINT FALLS IN THE CELL 
GRID CELL number: ROW 4» COLUMN 1C 0 <X< 240; 390 <Y< 650) 


AZIMUTHS 

CLASS LENGTH CIN mLE(S)» 

270 . 0 - 260.0 

C .94 

>xxxxxx 

260 . 0 - 290.0 

C .31 

> XX 

290 . 0 - 300.0 

U .96 

>xxxxxx 

300 . 0 - 310.0 

2.91 

>xxxxxxxxxxxxxxxxxxx 

31040 - 320.0 

2.59 

>xxxxxxxxxxxxxxxxx 

320 . 0 - 330.0 

0.89 

>xxxxx 

330 . 0 - 340.0 

1 .48 

>xxxxxxxxx 

340 . 0 - 350.0 

1.27 

>xxxxxxxx 

350 . 0 - 360.0 

0.52 

>xxx 

0 . 0 - 10.0 

1 .50 

>xxxxxxxxxx 

10 . 0 - 20*0 

1 .01 

>xxxxxx 

20 . 0 - 3 C .0 

2.46 

>xxxxxxxxxxxxxxxx 

30 . 0 - 40.0 

3.64 

>xxxxxxxxxxxxxxxxxxxxxxxxx 

40 . 0 - SO.O 

2.32 

>xxxxxxxxxxxxxxx 

SO.O- 60.0 

1 .98 

>KXXXXXXXXXXXX 

60 . 0 - 70.0 

0.69 

>xxxx 

70 . 0 - dO.O 

0.39 

> XX 

80 . 0 - 90.0 

0.47 

>xxx 

TOTALS 

26.54 

EACH X = 0.15 MILE 


RANDOMLY OIST'TiaUTED DENSITY DATA PROS- = 0.0 


NUMERICAL FREQUENCY 


2 >♦* 

1 >♦ 

4 >*♦*♦ 

7 >*«*«««« 

9 

3 >♦♦♦ 

5 >«««Hc4c 

3 >♦♦♦ 

2 >♦♦ 

4 >*♦*♦ 

4 >« 44 * 

9 

11 

9 

3 >♦*♦ 

1 >♦ 

2 >♦* 


BS EACH * ^ 1 UNITS 

RANDOMLY DISTRIBUTED FREO. DATA PRDB. - 0.7a56e^02 


Figure 1. Example of line printer output from the AZMAP program. Frequency and density 
histogram and Chi Square test options have been exercised. Grid cell coordinates 
are given in- millimeters and size is given, both in millimeters and the actual map 
units chosen. 



ROSE Program 


ROSE uses standard CalComp subroutines and hardware, producing rose 
diagrams from data generated on punched cards or magnetic tape in 
AZMAP. Those installations lacking this type of plotter or using other 
software subroutines will either have to forgo this last program, relying on 
the frequency-azimuth histograms, or restructure the program for use with 
their particular system. The program will not accommodate three-dimensional 
data such as jointing. Other programs are in the literature for stereonet 
plots^» Figure 2 represents a plot generated by this program. Less than 
40 seconds were required for processing 49 rose diagrams on an IBM 360/67 
computer with an additional 5 minutes for plotting on a CalComp 780 30-inch 
Drum Plotter. A source listing, abstract, and detailed instructions for use 
are included in the Appendix. 

ANALYSIS APPLICATIONS 

Discussion 

The following will illustrate several applications of these programs to a 
fracture analysis problem, A fracture trace analysis was performed on a 
225 square mile portion of northwestern Nolan and southwestern Fisher Coun- 
ties, Texas. Figure 3 is a geologic map of the area, located on the eastern 
shelf of the Permian Basin, It consists of a low relief (<200 feet) terrain 
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Figure 2. Cal Comp-generated rose diagram plot of 
fracture trace patterns. 



TEXAS STUDY AREA 


100 ® 35 ' 100 ® 30 ' 




Figure 3. Geologic map of the Texas study area^®> 

T - Tertiary Ogallala silts, sands, and gravel 
overlying predominately Cretaceous 
limestones of the Trinity and Fredricksburg 
Groups 

K - Cretaceous limestones of the Trinity and 
Fredricksburg Groups 
P-Tr - Permo-Triassic redbeds of gravels, 
sandstones, siltstones, and shales. Some 
gypsum horizons. 
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containing flatlying consolidated (Permian through Cretaceous) and unconsoli- 
dated Tertiary sedimentary materials. The fracture traces were mapped on 
1:20000 scale aerial photographs in flightlines and then transferred to 1:24000 
scale U, S, Geological Survey topographic maps by standard photogrammetric 
techniques. Figure 4 shows the fracture trace map. The number of observa- 
tions for this small area illustrates the need for automated data processing. 

The data were then digitized as discussed earlier and processed by the 
TRANSFORM program. Output consisting of 1486 fracture traces was 
subsequently used in AZMAP. 

Repetitive manipulations of the fracture data using the various algorithms 
within AZMAP, while varying the size of the grid cell, allowed an estimate 
of the grid cell size which provided the most efficient, consistent, and sufficient 
grid cell size and summarization technique. Based on these findings, a grid 
cell size of approximately 3. 25 x 3. 75 miles area was used. To gain further 
information about the variability of the fracture pattern, summarization was 
made at 1/2 cell increments. The following is excerpted from Podwysocki^^* 
The reader is referred to these publications for more detailed treatment of 
the data. 

Trend surface analysis^*^ of the number of fracture traces per grid cell 
summarized by AZMAP indicates that the third order surface shows fewer 
fracture traces over the Cretaceous-Tertiary rocks whereas greater 
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FRACTURE TRACE MAP NW NOLAN & SW FISHER COUNTIES, TEXAS. 


100 ° 40 ' 100 ° 35 ' 100 ° 30 ' 100 ° 25 ' 

■ I I I 



Figure 4. Fracture trace map of portions of Nolan and Fisher 

Counties, Texas. The upper and left axes correspond 
to the top and left margins of the rose diagram plot 
in Figure 2 and maps of Figures 3 and 5. 
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frequencies are encountered over the Permo -Trias sic clastic rocks. This 
may be interpreted as indicating 1) the older rocks are more fractured, 

2) the thin Tertiary mantle (usually less than six feet) has acted as a masking 
effect, 3) the rook types respond differently to differences in mechanical 
properties, or 4) two stress fields may have operated on the older rocks 
while only one affected the younger. Higher-order trends tend to align 
themselves with the direction of flightlines, a bias most likely attrib\itable 
to changes in operator accuity. This suggests tiiat photographs should not be 
mapped in flightlines, but that they should be randomly selected to try to 
eliminate or distribute this bias evenly during the project. 

Cluster analysis, a multivariate analysis technique used to classify samples 

15 16 17 

on the basis of their similarity or variability, is becoming widely used * * 

A cluster analysis program using covariance measures written by Rubin 
18 

and Friedman was applied to the frequency- azimuth histograms generated 
by AZMAP in order to classify areas of like-behaving fracture patterns. 

Each grid cell was treated as a sample consisting of 18 variables (18 azimuth 
classes of 10 degrees each). Optimum grouping was determined by plotting 
the program-generated optimization measure (log|T| /|w| ) against the 
number of groups, which was varied from two through twelve. The greatest 
rate of increase in the measure was noted between the two and three group 
levels, with additional minor inflections occurring at higher grouping levels. 
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Based on the three-group classification, the bedrock geology of the area could 
be mapped (Figure 5). 


TEXAS STUDY AREA 
100 ® 35 ’ 100 ® 30 ' 



Figure 5. Results of a three groi:^} cluster analysis classification 
of rose diagram patterns observed in Figure 2 based on 
covariance measures. Group 1 - mainly Cretaceous 
and Tertiary rocks; group 2 - transition between groups 1 
and 3; group 3 -mainly Permian and Triassic rocks. 
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Conclusions 


These programs are an attempt to facilitate easier handling of fracture trace 
or lineament data in order to provide a basis for the application of rigorous 
statistical treatment. Users with a knowledge of FORTRAN can, in addition, 
modify the programs to suit their individual needs. The rose diagram program 
also facilitates the presentation of the fracture data in a suitable form for 
qualitative interpretation of the spatial relationships between areas for tectonic 
analysis. Additional numerical techniques applicable to fracture analysis have 
been suggested. 
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APPENDIX 


PROGRAM SOURCE LISTINGS 



VECTOR TRANSFORM PROGRAM 


Ptece4ii>8 

A-3 



c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


TRFMOOlD 

VECTOR T^(kNSFORM PROGRAM TPFM0020 

TRFM00 30 
TRFM0040 

THE PROGRAM WAS WRITTEN BV MELVIN PODWYSQCKt OF THE GEOSC 1 ENCESTPFM0050 
OEPTo, THE PENNSYLVANIA STATE UNIVERSITY, APRIL, I R72 FOR USE TRFM0060 
ON THE T0M 360/67 COW^UTER, AND WAS MODIFIED TN ARRILa 1973 TRFM0070 
FOR use ON OTHER COMPUTERS HAVING THE EQUIVALENT CORE STORAGE TRFMOOSO 
OF 78K BYTESa TRFM0090 

TRFMOlOO 

PROGRAM CALCULATES VECT39 LENGTH (VECLEN), AZIMUTH (VECAZMI, TRFMOllO 

SLOPE <AJ, y INTERCEPT (BJ AND X & Y MIDPOINTS (XMID & YMfD> FOR TRFM0120 
EACH VECTOR GIVEN ITS 3EGINNING IXl.YM AND END { X2 , Y2 ) POINTS TRFM0130 
IN MAP SPACE* PROGRAM ASSUMES X IS ♦ TOWARDS THE RIGHT AND Y TRFM0I40 

IS + DOWNWARD. DUE NORTH“SOUTH AND EAST-WEST DATA ARE TREATED AS TRFM0150 
SPECIAL CASES. CONTROL CAROS ARE READ FROM THE CARD READER WHILE TRFM0160 
DATA CAPOS ARE READ FROM ANY UNIT DECLARED BY »XTAPEl* ON CONTROLTRFMO 1 70 
CARO 3. DATA CARDS ARE GENERATED ON ANY UNIT DECLARED BY • I T APE 2« TRFMO 1 80 
IN A format ACCEPTABLE TO "A2MAP** PROGRAM, COORDINATES MUST BE TRFMD190 
IN millimeters, if NOT, CONTROL CARD 3 ALLOWS A TRANSFORMATION TOTRFM0200 
BE PERFORMED. PRINTED OUTPUT IS ALSO PRODUCED. TRFM0210 

TRFM0220 

(0,0) X ♦ > TRFM0230 

TRFM0240 

TRFM0250 
TRFM0260 
TRFM0270 
TRFMO20O 
TRFM0290 
TRFM0300 
TRFM0310 
TRFM0320 

I I TRFM0330 

V I TRFM0340 

TRFM0350 
TRFM0360 
TRFM0370 

ALL NUMERIC INPUT DATA IS RIGHT JUSTIFIED; INDICATES INTEGER TRFM0380 

FORMAT, "F” INDICATES Floating point format, ••A" indicates CHARAC-TRFM0 390 


ALL VECTORS 
LIE IN THIS 
AREA 


NORTH 

A 

I 

I 


TER FORMAT, PRECEEDiNG NUMBERS INDICATES COLUMNS USED FOR 

EACH PARAMETER. TO SPECIFY NONUSE OF AN 3PTION, PUNCH 0 OR 
LEAVE BLANK. 

A«r*^#?*i>ACONTROL CARD 1 --— — TITLE CARD 


TRFM0400 
TPFM0410 
TRFM0420 
TRFM0430 
TRFM0440 

TI1LE WILL BE PRINTED AT BEGINNING OF PRINTED OUTPUT (20A4, 41-60) TRFM0450 
♦ ^ 1 44# #4^ CONTROL CARD 2-- — V AR I ABLE INPUT FORMAT CARD TRFM0460 

TO READ X 8 Y VALUES FROM CARDS. FORMAT MUST BE ENCLOSED IN TRFM0470 

PARANTHESES and BEGIN IN 41. VALUES SHOULD BE READ IN THE FOLLOW- TRFM0480 
ING order; XL,Y1,X2,Y2 C20A4,4l-80> TRFM0490 

*444444*fC0NTP0L CARO 3 —OPTIONS CARD TRFMOSOO 

ITAPEl=LOGTC AL UNIT FOR READING DATA CAROS (I2,41<-2) TRFM0510 

IT APE2=L0GICAL UNIT FOR WRITING OUTPUT USEABLE BY ••AZMAP'* <12, TRFM0520 
4 3-41 TRFM0530 

NPRINT--PRINT COMMAND F3 ^ OUTPUT (11,1 IN 45) TRFM0540 

NPUNCH — PUNCH COMMAND FOR OUTPUT ACCORDING TO FORMAT FOR "AZMAP« TRFM0550 

PROGRAM ON ANY UNIT DECLARED BY ’ITAPEa* (11,1 IN 46) TRFM0560 

NTRAN TRANSFORMATION COMMAND TO CONVERT UNITS DF INPUT DATA TO TRFM0570 

MILLIMETERS (II. 1 IN 47) TRFM0580 

note; THE FOLLOWING PARAMETER IS USED ONLY IF *NTRAN* IS PUNCHED TRFM0590 
1. UNITS WILL BE CONVERTED TO MILLIMETERS FOR USE IN ••AZMAP”. TRFM0600 














A -5 









uuu ouu uuu uuuuuu 


C CONVERTED COORDINATE. VALUE = ORIGINAL VALUE » *CONV • 

C CONV=TPANSFOPMATION FACTOR TO CONVERT UNITS ON INPUT DATA CAROS 

C TO MILLIMETERS <F6.2,#B-13l 

C CAROS— 

C READ FROM ANY UNIT DECLARED BY MTAPE1«« END OF DATA CARDS IS 

C SIGNIFIED WHEN THE PROGRAM ENCOUNTERS A CAPO WHOSE VALUES OF 

C Xl*Yl,X2 e Y2 APE ALL 0«0 

C 

DIMENSION FMTRO< 20>*TIT.E<20) 

DATA IREAD/S/. IPRINT/6/ 

C 

C READ CONTROL CAROS 

C 

REA0<IREAD»5) < T ITLE< L ) . U » 1 * 20 ) 

READ Cl REA0*5 ) ( FMTROCL ) • L = 1 • 20 > 

REAOCIREAO*? ) I TAPE 1. ITAPE2* NPR INT. NPUNCH#NTRAN *COMV 
WRITEC IPRTNT,^ > CTITLECL > *L=1 *20) • (FMTROCL » #L=1 ,20) 
t FCNTRAN* E0» 1 ) WR ITE < IPR I NT * 1 0 > CONV 
TFCNPRINT.EQ# 1 ) WRITE<l=»R INT* U) 

READ X t Y COORDINATES FROM UNIT MTAPEl* 

NUMjeO 

12 PEAOCITAPEl.FMTRO ) Kl,YltX2«Y2 

IFCXl* E0«0**ANDa Yl«EQ«Oo » ANO« X2oEQ«0»*ANO*Y2«EO* Oo » GO TO 100 
NUM=£NUM-M 

CONVERT DATA TO MILLIMETERS IF NECESSARY 

IFCNTRAN) 20 t20» 15 
15 Xl-XIACONV 
Yl5=Y lit CONV 
X2*X2* CONV 
Y2=Y 2* CONV 

ORDER BEGINNING AND END COORDINATES FOR "AZMAP« PROGRAM 

20 IF(Y 1« EO» Y2« ANDo X1.LT*X2I GO TO 22 
IFCY1-Y2) 21*23«23 

21 22«Y2 
Y2=Y1 
Y1=Z2 

22 Zl=X2 
X2=X1 
XI =2 1 

CALCULATE VECLEN 

23 VECLEN-SQRTC < ( X2- X I Mi*2 ) 4^ C ( Y2-Y1 
CALCULATE SLOPE <A* 

IFCX2-XI) 24t25« 24 

24 A = CY2-Yl )/(X2-XX J 
GO TO 30 

25 A^-573o0 
VECA ZM = 0« 0 
0^ 50 0 0 00* 

GO TO 50 

CALCULATE VECAZM 


TRFM0610 

TRFM0620 

TRFM063O 

TRFM0640 

TRFM0650 

TRFM0660 

TRFM0670 

TRFM06SO 

TRFH0690 

TRFM0700 

TRFM0710 

TPFM0720 

TRFM0730 

TRFM0740 

TRFM0750 

TRFM0760 

TRFM0770 

TRFM0T90 

TRFM0790 

TRFMO0OO 

TRFMOSIO 

TRFM0920 

TRFM0B30 

TRFM0540 

TRFMO05O 

TRFM0850 

TRFM0870 

TRFMoeeo 
TRFM0B90 
TRFM0900 
TRFM0910 
TRFM0920 
TPFM0930 
TRFM0940 
TRFM0950 
TRFM0960 
TRFM0970 
TRFMO90O 
TRFM0990 
TRFMIOOO 
TRFMIOIO 
TRFM1020 
TRFM1030 
TRFH1040 
TRFM1050 
TRFMl 060 
TRFM1070 
TRFMl 060 
TRFM1090 
TRFMl 100 
TRFMl 110 
TRFMl 120 
TRFMl 130 
TRFMl 140 
TRFMl I 50 
TRFMl 160 
TRFMl 170 
TRFMl 180 
TRFMl 190 
TRFM1200 
TRFM1210 
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30 1F(A> 35.33#33 

33 y/ECAZM = 2T0®0 ^ < ( ATAN( A > • 1 80<» ) / 3« 1 4 1 59 ) 

GO TO 40 

35 \/ECAZW = 90« -» ( ( AT AN ( ABS< A M4 1 80» I y 3« 1 4 159 ) 

CALCULATE Y INTERCEPT •B* 

40 B={( C Y X2 1/ ( X2- XI n ^Y2 

CALCULATE MIDPOINT OF VECTOR 

so xMiD=( xa^xi I /a«o 

ymio=( V2+vn/2*0 


OUTPUT 

I F <NPR 
62 WPITE( 
75 IF(NPU 
77 WRITE( 
GO TO 
100 WRITE< 
5 FORM AT 
7 FORMAT 
9 FORMAT 

10 FORMAT 

11 FORMAT 

1 C A 7M ♦ o 

65 FORMAT 
AO FORMAT 
105 FORMAT 
STOP 
END 


INT) 76. 75,62 

TPRINT.eSl XI , YJ ,X2. Y2.VECLEN.VECAZM* A.B.XMID.VMTO 
NCH) 12.12.77 

ITAPE2. eOl X1.Y1.X2. Y2. VECLEN . VEC A Z M . A .B . X MIO . YMT O 
12 

IPRINT.1051 NUM 
(20A4 ) 

(212. 31 1«F6«2) 

C IHI, 20A4/. IH 0. 35HV ARl ABLE INPUT FORMAT FDR X t Y IS 
(1H0.44HALL VALUES WILL BE MULTIPLIED BY A FACTOR OF.F 
( IHO. 4X.2HXI . IIX. 2HYI . I 1 X.2MX2. 1 1 X. 2HV2 ,9 X. 6HVECLEN.5 
XtlHA. 13X. IHB* 12X.4HXMID. 9X. 4HYM1D/,^| 

UH .5CF7, 1.6X ).F5.1 • 6X . F 6* 3 • 6X «F9* 1 • 2 ( 6X, F7* 1 1 ) 

( 6F6* l.F9*4. F10> 1 . 2F7* 1 ) 

( IHO. 20HNUMBER OF VECTORS « .1101 


TRFMl 220 
TRFM1230 
TRFMl 240 
TRFMl 250 
TRFMl 260 
TRFM1270 
TRFMl 2B0 
TRFM1290 
TPFM1300 
TRFM1310 
TRFMl 320 
TRFMl 330 
TRFMl 340 
TRFMl 350 
TRFM1360 
TRFMl 370 
TRFMl 3«0 
TRFM1390 
TRFM1400 
TRFMl 410 
TRFMl 420 
TRFMl 430 
TPFM1440 
TRFM1450 
TRFM1460 
« 20 A41TRFMI 470 
12*5) TPFM140O 
X .6HVETRFMI 490 
TRFMl 500 
TRFM1510 
TRFM1620 
TPFM1530 
TRFMl 540 
TRFM15S0 
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AZMPOl 00 
AZMPO 110 
AZMP0120 
A2MP0 1 30 
AZMPO 140 
AZMPOl 50 
AZMPOl 60 


C file A AZMPOOIC 

C A7MAP PROGRAM AZMP0020 

C ♦♦****ritA**4»**»#i*^<ii»!ik*4*^^^*** AZMP0030 

C AZMP0040 

C THE PROGRAM WAS WRITTEN BY MELVIN PODWYSOCKT OF THE GE OSC I ENCES AZMPO 0 60 

C DEPT*. THE PENNSYLVANIA STATE UNIVERSITY « APRIL# 1072 FOR THE AZMP0060 

C IBM 360/67 COMPUTER. AND WAS MflDIFIEO IN APRIL. 1973, FQP USE AZMP0070 

C ON OTHER COMPUTERS HAVING THE EQUIVALENT OF 160K BYTES ST OP AG E# AZMPO 0 00 

C AZMP0090 

C PROGRAM SUMMARIZES FREQUENCY O I STR I BUT 1 0N5 OF VECTOR DATA IN 

C VARIABLE MAP GRID AND AZIMUTH CLASS SIZES. PROGRAM ALLOWS JP TO 

C 90 AZIMUTH CLASSES. FROM 270 THRU 0 TO 90 DEGREES. SUMMARIZING 

C DATA AS TOTAL LENGTH OF V ECTORS /A ZI MUTH CLASS (DENSITY) OR NUM- 

C BER OF VECTOPS/AZIMUTH CLASS (FREQUENCY). UP TO 2000 VECTORS 

C MAY BE USED. X AXIS IS 4 TO RIGHT t Y AXIS IS 4 DOWNWARD. DATA 

C ARE READ FROM CAROS GENERATED BY VECTOR TRANSFORM PROGRAM. 

C CONTROL 6 TITLE CAPOS A^E READ FROM CAPO READER WHILE DATA CARDS AZMP0170 

C MAY BE READ FROM ANY UNIT DECLARED BY MTAPE2* ON CONTROL CARO I.AZMP0160 

C AZMP0190 

C ALL NUMERIC INOUT DATA IS RIGHT JUSTIFIED; "I** INDICATES INTEGER AZMP0200 

C FORMAT. ”F»* INDICATES FLOATING POINT FORMAT. **A •• INDICATES CHA AZMP0210 

C RACIER FORMAT. PRECEEOING NUMBERS INDICATES COLUMNS USED FOR AZMP0220 

C EACH PARAMETER. TO SPECIFY NONUSE OF AN OPTION. PUNCH 0 OP LEAVE AZMP023O 

C BLANK# AZMPO 240 

C AZMP0250 

C *4 |i**> CONTROL CARO OPTIONS CARD AZMP0260 

C X INC = INCREMENT OF X-AXIS TRAVERSE IN MILLIMETERS (14. 41-4) AZMP027O 

C Y 1NC = INCREMENT OF Y-AXlS TRAVERSE IN MILLIMETERS (I4.45-8) AZMPO20O 

C XSTART = STARTING POINT FD ^ X-AXIS TRAVERSE IN MILL I METER S ( I 4 . 9- I 2 > AZMP0290 

C YSTARTssST APT ING POINT Y-AXlS TRAVERSE IN Ml L LI ME TERS ( I 4 . 41 3- 16 AZ MPO 300 

C XSTOPsENO of X-AXlS TRAVERSE IN MILLIMETERS (14.417-20) AZMP0310 

C YSTOP=END OF V-AXIS TRAVERSE IN MILLIMETERS (14.421-24) AZMP0320 

C note: program successively scans DATA IN MAP GRID CELLS • XCELL ♦ AZMPO 3 30 

C BY 'VCELL* IN SIZE. INCREMENTING BY •XINC* UNTIL *KMAX* > AZMP0340 

C *XST0P*. when •riNC* is incremented. program terminates when AZMP0350 

C 'YMAX* > *YSTDP*. NONE OF THE ABOVE 6 VALUES CAN BE NEGATIVE. AZMP0360 

C AZCLAS=^AZIMUTH CLASS WIDTH IN DEGREES ( F 3. 1 . 425-27) AZMP0370 

C AMPSCL*MAP$C ALE IN UNtTS/SCALE; SEE SCALE PAPAM» BELOW (F5.4.426- AZMPO30O 

C 32) AZMP0390 

C KTYPE--SEL6CTS SUBPOUTIN- FOR CLASSIFYING VECTOR DATA <12.433*34) AZMP0400 

C PUNCH 1 FOR SUBR. PAt>T. CONSIDERS ONLY THAT PART OF VECTOR AZMP04 10 

C WHICH LIES WITHIN THAT CELL AZMP0420 

C PUNCH -1 FOR SUBR. MID. CONSIDERS WHOLE VECTOR IN CELL IF ITS AZMP0430 

C MIDPOINT FALLS WITHIN THAT CELL AZMP0440 

C NOHTSr PRINT COMMAND FO^ DENSITY HISTOGRAM (Tl.l IN 435) AZMP0450 

C NFMIST — PRINT COMMAND FOR FREQUENCY HISTOGRAM (Tl.l IN 436) AZMP0460 

C NPUNCH — PUNCHED CARD OUTPUT FOR EACH GRID CELL AZIMUTH DISTRIBU- AZMP0470 

C r ION IS GENERATED AS >ER CONTROL CARD 2A BELOW (II. I IN 437) AZHPQ40O 

C 0HINC=NUMERICAL VALUE OF EACH 'X* INCREMENT (ACCORDING TO SCALE AZMP0490 

C PARAMETER; SEE BELOW) FOR DENSITY HISTOGRAM ( F5.2 « 43 8-42 ) AZMPQ500 

C NFHl NC=NUMER ICAL VALUE OF EACH INCREMENT FOR FREQUENCY HISTO- AZMPOSIO 

C GRAM ( 13. 443-45) AZMP0520 

c SCALE5:SCALE units for PRINTOUT ( I • E« MI lE S , M M. .ETC. ) . ( 2A4 . 446 - 53 ) AZMP0530 

c note: when vectors are measured on a 1:24000 SCALE MAP AND 0UT-AZMP0540 

C PUT IS DESIRED IN MILES . • AMPSCL • = • 01 49 ( I . E. 1 MM.=.0149 MILES I A Z MP0550 

C I T APEl^LOGICAL UNIT FOP READING DATA CARDS GENERATED BY «TRANS“ AZMP0560 

C FOPM« PROGRAM (12.454*55) AZMP0570 

C NCMI — COMMAND TO TEST FREQUENCY AND DENSITY DATA FOR RANDOMNESS AZMP05BO 

C IN EACH GRID CELL (II. I IN 460) AZMP0590 

C XCELL=CELL SIZE IN X DIRECTION (T4. 461-64) AZMP0600 


PRECEDING PAGE BLANK NOT FILI/IED 


A-ll 



nnn nnrv nor* 


YCELL=ceLL SIZE IN Y DIRECTION (14,4^65*60) AZMP0610 

NUN=tNUHBER OF VECTORS (AS GIVEN IN TRANSFORM PROGRAM) (I4*«S9*7ZI AZMP0620 
A>***»*AC0NTP0L CARO 2--— --TITLE CARO AZMP0630 

TITLE WILL 0E PRINTED AT TOP OF EACH HISTOGRAM C20A4«Al-80) A7MP0640 

4i*:*AA*** CONTROL CARO 2A-^- - * V AR I A0LE OUTPUT FORMAT CARO AZMP0650 

TO BE USED ONLY IF NPUNC-I (CONTROL CARO 1) IS PUNCHED 1 AZMP0660 

PUNCHED OUTPUT WILL CONSIST OF A CELL ROW. COLUMN. X S Y MI DPOl NTS AZMP0670 
OF CELL. AND THE VALUE OF EACH AZIMUTH CLASS SUMMATION AS PER AZMP0680 
VARIABLE "ITYPE*. THE FOLLOWING FORMAT IS SUGGESTED! 

( 21 4 . 2F 7« 2. 0F 7. 2/1 OF7. 2 ) 

NOTES THE LAST SFT DF VARIABLES MUST CORRESPOND TO THE NUMBER 
OF AZIMUTH CLASSES <I*E* I 8 0 • / • A ZCLAS * = NUMBER OF CLASSES)® 

IN THE ABOVE EXAMPLE IT-S 18. FOR DENSITY DATA IT SHOULD BE AS 
ABOVE. BUT IN INTEGER FORMAT^ loE* < • •• • • • .6 15/10 15 ) 

FMTPCH-*OUTPUT FORMAT AZIMUTH CLASS DATA® MUST BE ENCLOSED 

IN PARANTHESES AND START IN Wl® < 16A4 * »l-72) 

ITYPE--PUNCH 2 FOR DENSITY DATA (TOTAL LENG TH/AZ IMUTH CLASS) 

PUNCH 1 FOP FPEQUENCY DATA (NUMBER OF VH CTQRS /A ZI MUTH CLASS) 

OUTPUT WILL BE IN UNITS SPECIFIED BY “AMPSCL* AND *SCALe» 

( II . A77) 

Ie^*^*OATA CAROS-- — * 

VECTOR DATA INPUT FROM VECTOR TRANSFORM PROGRAM 


4*4 


AZMP0690 
AZMP0700 
AZMP0710 
AZMP0720 
AZMP0730 
AZMP0T40 
AZMP0750 
AZMP0760 
AZMP0770 
AZMP0760 
AZMP0T90 
AZMPOeOO 
AZMPOetO 
AZMP0820 
AZMP0830 
AZMPQ840 

DIMENSION TITLEC 20). SCALE ( 2).FMTPCH( 16) . 23( 90) AZMP08S0 

COMMON XI (2000). Yl(2 000) «X2(2000 ). Y2( 200 0 ) .VECAZMf 2000) . VECL EN ( 2 00 AZ MPO 8 60 
10 ) »A ( 2000 ).B (2000). XMIDC 2 000) .YM ID ( 2000) .CLAMlN( 90 ) .CLAMAXCROI . AZLAZMP0 870 
2EN(90) .NAZFRQC90I tXMTN.YMIN.XMAX.YMAX.FLAG. VLEN. I HUM AZMPO08O 

INTEGER XST3P* YSTOP, XiNC. YINC. X STAR T,YST ART ,XMlN ,Y «I N . X MA X , Y M AX . FL AZ MPO 890 


lAG.XCELL. YCELL 

DATA I XS/IHX/. ISTARS/1K4/. lREAD/5/. I PPl NT/6/ . 1 PUNCH/7/ 

READ CONTROL INPUT INFORMATION AND TITLE CARO 

READ (I REA0.5 ) XI NC. Y INC. XSTART. YSTAPT. XSTOP .YSTOP. AZCLAS . AMPSCL ♦ 
XXTYPE.NDHIST.NFHIST.NPUMCH.OHINC.NFHINC. < 5CALE(L|.L=1 .2). ITAPEl. 
2NCHI . XCELL. YCELL .NUM 
READ ( I READ. 6 ) ( T ITLE (L ) .L- 1 . 20 ) 

IF(NPUNCH.EQ® 1 I R E AO ( IRE AO . 7 ) (FMTPCH(L) «L=1 .18) . I TYPE 
Z1-XCELL4 AMPSCL 
Z2=Y CELL* AMPSCL 
IF(KTYPE) 19.9.19 
9 WP.ITE( IPRINT. 10) 

STOP 

READ DATA GENERATED BY VECTOR TRANSFORM PROGRAM 
19 DO 25 1^1 .NUM 

READ (I TAP El. 20 ) X 1 ( T I . YU I > . X2 ( I ). Y2( I t. VECLENC I ) . VECAZMC I ) . AC I I , 
1 8( I) «XM IO( I) ,YMID( 1 I 

25 CONTINUE 


GENERATE AZIMUTH CLASSES 


AZMP0900 
AZMP091 0 
AZMP0920 
AZMP0930 
AZMP0940 
AZMP0950 
AZMP0960 
AZMP0970 
AZMP0980 
AZMP0990 
AZMPl 000 
AZMOl 010 
AZMPI020 
AZMP1030 
AZMP1040 
AZMPl 050 
AZMP1060 
AZMPl 070 
AZMPIOBO 
AZMPl 090 
AZMPl ICO 
AZMPl I 10 
AZMPl 120 
AZMPl 130 
AZMPl 140 


30 CLAM IN ( I ) =270o0 
1=1 

NCLASS=l 

CLAMAXC n=CLAMlN(l )♦ AZCLAS 

40 IF(CLAMAX( 1 )®GT:, 270, O. AND. CLAMAXC T >®LE®360. 0 ) GO TO 41 
TF<CLAMAX( I )-90,0) 41.60.60 

41 I = U^1 


AZMPl ISO 
AZMPl 160 
AZMPl 170 
AZMPl 180 
AZMPl 190 
AZMPl 200 
AZMP1210 
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NCLASS-I 

CLAM IN { I >5=CL AMAX ( I- I » 

CLAM AX C I ) s= CL AH INI I 1 + A2CLAS 
IF(CLAMAX< n-360« 01 40*50iS0 

50 NCLASS = I*‘l 
T=NCLASS 

CLAM iNf I>=CLAMAXC I-l )- 360<, 0 
CLAM AX ( I I =CL AM IN C I J ♦ AZCLAS 
GD TO 40 
C 

C SCAN & SUMMARIZE DATA FD^ EACH GRID CELL 

C 

60 00 300 YM IN* YSTARTf YSTOR* Y INC 
YMAX=YMIN+YCELL 
IFCYMAX-VSTOP) 62»62t400 
62 DO 300 KM !N= XSTART* XSTOR# XINC 
XMAX=XMIN4XCELL 
I F<XMAX-^XSTQP I 64»64.300 

64 00 65 L*1 .90 
AZLEN(L) = 0^ 

65 NAZFRQ<LI=0 
WVLEN*0. 

NEWFRQsO 

DO 170 1=1. NUM 

INUM=I 

FLAG=0 

IFCKTYPEI 6fl. 9, 66 

66 CALL PART 
GO TO 69 

69 CALL MID 

69 IF(FLAG>70.75.70 

TEST FOR E • M DATA 
C 

70 IF ( VECAZMCI )-270« ) 76*71.76 

71 1FCFLAG)72.170.73 

72 WVLEN=WVLEN+VECL£N<I } 

GO TO 74 

73 WVLEN=WVLEN+VLEN 

74 NEWFRQ=NEMFRQ^1 
C 

C ADO RESULTS TO APPROPRIATE AZIMUTH CLASS 

C 

75 IFU-NUM) 170*76.7© 

76 IFCI-NUMl 62.78.70 

78 IFCAZLENC 1 )- AZLEN {NCLASS I ) 79 * 80.80 

79 A ZLEN< NCLASS 1 = AZLEN< NCLASS I +WVLEN 
NAZFRQINCLASS > =N AZFR 0 C NCL A SS > ♦NEWFRQ 
GO TO 81 

80 AZLENC 1 )= AZLEN( 1 M-WVLEN 
NAZFRQt n=NAZFRQ < 1 ) + NEITFRQ 

91 IF<VeCAZM(n«NE*2T0««AND*FLAGoNE.0oAN0*T*EQ*NUM| GO TO ©2 
GO TO 170 

82 DO 85 J=l. NCLASS 

IF ( VECAZHCI |« GF«CLAM1N( J ) «AND« VECAZMin o LT« CLAM AX< J M GO TD ©4 
GO TO 05 
94 NTVPE=J 
GO TO 90 
85 CONTINUE 

90 IFIFLAG) 100* 170. 150 
100 AZLEN( NTYPEI=AZLEN(NTYPEM-VECLEN<I > 


AZMP1220 
AZMPl 230 
AZMP1240 
AZMPl 2 50 
AZMP1260 
AZMPl 270 
AZMPl 280 
AZMP1290 
AZMPl 300 
AZMPl 310 
AZMP1320 
AZMPl 330 
AZMPl 340 
AZMPl 350 
AZMPl 360 
AZMP1370 
AZMPl 380 
AZMPl 390 
AZMP1400 
AZMP1410 
AZMP1420 
AZMP1430 
AZMPl 440 
AZMP1450 
AZMP1460 
AZMP1470 
AZMP14B0 
AZMPl 490 
AZMP1500 
AZMP1510 
AZMPl 520 
AZMPl 530 
AZMPl 540 
AZMP1550 
AZMPl 560 
AZMP1570 
AZMP1580 
AZMPl 590 
AZMPl 600 
AZMP1610 
AZMPl 620 
AZMP1630 
AZMP1640 
AZMPl 650 
AZMPl 660 
AZMPl 670 
AZMP1680 
AZMP1690 
AZMP1700 
AZMP1710 
AZMPl 720 
AZMPl 7 30 
AZMP1740 
AZMPl 750 
AZMPl 760 
AZMPl 770 
AZMP1780 
AZMPl 790 
AZMPl 800 
AZMPieiO 
AZMPl 820 
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150 
1 60 
170 


1 ao 

lai 


152 


154 


195 

196 
195 

214 

215 

217 
21 R 

219 

220 


GO TO 160 

AZLeN< NTYPEI = AZLEN{NTYPE> <-VUEN 
N AZFRQt NTYPE I = NA 2FR0 ( NTYPE ) + l 
CONT INUE 

OUTPUT 

TOTLN^O# 

NFRQ-0 
NXERR=0 
MAST EP=0 

DO 180 N=l,NCLASS 
T OTL N=^T OT LNf A ZLE N < N ) 

NFRQ=NFRQ*NA ZFRCMNl 
CONT INUE 

TEST FOR RANDOMNESS OF AZIMUTH D ISTR I BUT I ON S/CELL BY CHl SQUARE 

IFCNCHI ) 184,184«iei 

CLASS^NCLASS 

FRQsJNFPQ 

OENE XP=TOTLN/CLA SS 
FRQE XP=FRQ/CL ASS 
OCS^sQ. 

FCS-0. 

DO 102 LCS=1*NCLASS 

DCS=0CS-M C AZUEN(LCS)-0ENEXP IA42I /DENEXP 

FCS=FCS-M CNAZFRQ (LCS 1-F9QEXP 2 )/FRQE XP 

CONT INUE 

NDFaNCLASS-1 

0CHPR6=PR BCHK OCSpNOF » 

FCHPRB=PR8CHI(FCS«N0F> 

PRINT DATA FOR EACH CELL 

HROtf YMIN4Y INC I /VINC 
NCOL=( XMIN+XINC) /XINC 
GXMI Ds < XM !N4>XMAX )/2^ 

GY MI D= ( YMIN4 YMAX 1/2* 

WRITE< I PRINT. 185 ) (TITLE<L ) .L=l.20 } 

WRITE! 1PRINT.1061 XCELL. Z 1 . ( SC ALE( L > .L^1.2I . YCELL. Z2 .< SCALE! L > .L 

I .2) 

IF(KTYPE) 195.9. 196 
WRITE! 1 PRINT. I 80 > 

GO TO 190 
WRITE! IPRTNT# 187 ) 

WRITE! IPRINT. 190 I 
WRITE! IPR INT 9 200 > 

00 230 T-l.NCLASS 
Z3!I >=AZLEN! I )4AMPSCL 

WRITE! IPR INT .2051 CL AMIN! I 1. CLAMAX! I I .Z3( II .NAZFRO!! I 
IF<NDM1ST >220.22 0.214 
NUMX-Z3< I )/DHlNC 
IF (NUMX >217.216, 217 
WRITE! tPR INT. 2161 
GO TO 220 

1 F!NUMX-50 1219. 219.218 
NUMX =50 
NXERR=1 

WRITE! IPR INT. 266) ( I XS . I UK A= 1 « NUMX ) 

IFCNFH 1ST 1230.230.221 


NROW.NCOL. XMIN. XMAX.YMIN. YMAX 
! SCALE! L >.L=1 .2) 


AZMP1630 
AZMP1840 
AZMPl 050 
AZMP1860 
AZMPl 870 
AZMPl 880 
AZMPl B90 
AZMP1900 
A7MP19I0 
AZMP1920 
AZMPl 930 
AZMPl 940 
AZMP1950 
AZMP1960 
AZMPl 970 
AZMP1980 
AZMP1990 
AZMP2000 
AZMP2010 
AZMP2020 
AZMP2030 
AZMP2040 
AZMP2050 
AZMP2060 
AZMP2070 
AZMP2080 
AZMP2090 
AZMP2100 
AZMP21 10 
AZMP2120 
AZMP2130 
AZMP2140 
AZMP2150 
AZMP2160 
AZMP21T0 
AZMP2180 
AZMP2190 
AZMP2200 
AZMP2210 
AZMP2220 
:1 AZMP2230 
AZMP2240 
AZMP2250 
AZMP2260 
AZMP2270 
AZMP2280 
AZMP2290 
AZMP2300 
AZMP2310 
AZMP2320 
AZMP2330 
AZMP2340 
AZMP2350 
AZMP2360 
AZMP2370 
AZMP2380 
AZMP2390 
AZMP2400 
AZMP2410 
AZMP2420 
AZMP2430 
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c 

c 

c 


221 NUMAST=NA2FRQ< I I/NFHINC 
IFfNUMAST 1224, 22 2. 224 

222 KRITE( 1PRINT^223I 
GO TO 2 30 

224 IF fNUMAST*>43l226«226 «225 

225 NUMAST=43 
MAST ER^I 

226 WRITE! IPRINT.207 I < I STARS . LI0R ALs=l . NUMAST I 
230 CONTINUE 

PUNCH CARD OUTPUT AS PER CONTROL CARD 2A 

IF(NPUNCH) 2 52 *2 52. 2 34 
234 IFfITYPE-1) 330,238,236 


A;MP2440 

AZMP2450 

AZHP2460 

AZMP2470 

AZWP2480 

AZMP2490 

AZMP2500 

AZMP2510 

AZMP2520 

AZMP2530 

AZMP2540 

AZMP25S0 

AZHP2560 

AZMP2570 


236 WRITE! IPUNCH.FMTPCHI NPOW , NCOL , G XM I D, GYM I O, < Z3( < OA 1 « KOA^^l • NCLASS I AZMP2580 
GO TO 252 AZMP2590 

238 WRITE! IPUNCH.FMTPCH) NROW, NCOL ,GXMI D ,GYMlO, C NA2F RQ CKQAI ,KQA^1 , AZMP2600 

INCLASS) AZMP2610 

AZMP2620 

PRINT SUMMARY DATA FOR EACH CELL AZMP2630 

AZMP2640 

252 WRITE! IPRINT, 2401 AZMP2650 

TLENM-TOTLN4 AMPS CL AZMP2660 

WRITE! tPRiNT, 2501 TLENM,NFRQ AZMP2670 

IF!NOHTSTI270,270,260 AZMP2650 

260 WRITE! IPPINT, 265 1 OH INC, ! SCALE! L ) , L= 1 , 2 1 AZMP2690 

IFINXERR.EO. 1 ) WRI TE ! IPR INT • 290 > AZMP2700 

270 IF(NFMIST 1296, 296,200 AZMP2710 

280 WRITE! IPRINT ,2851 NFHINC AZMP2720 

IF(NASTER«E0«1 J WRITE! IPR I NT, 295 1 A2MP2730 

296 IFCNCHUEO*! I WR ITE ! IPR I NT , 297 ) OCHPRB,FCHPR0 AZMP2740 

300 CONTINUE AZMP2750 

5 FORM AT! 614, F3« 1,F5«4,I2,3I 1 ,F5* 2 , 1 3 ,2A4 , I 2, 4X ,1 1 ,3141 AZMP2760 

6 F0PMAT!20A4> AZMP2770 

^ FORMAT ( 18 A4,4X,I 1 1 AZMP2780 

10 FORM AT (lHl, 30 X,eiHNEITHER SUMMARIZATION TECHNIQUE {PART OR MIDI WAAZMP2790 


IS SPEC IFIEO»*»*** JOB AB0RTED44444I 
20 F0PMAT!6F6,1 ,F9« 4,FI 0*l,2F7«t I 
IBS FORM AT! 1H1,20X,20A4> 


AZHP2800 

AZMP2B10 

AZMP2820 


186 FORM AT ! IH0,20X, 18HEACH GRID CELL IS ,I4,6H MM# C , F7 • 3 , 2A4, SH 1 BY AZMP2830 
t,l4,6H MM« ( ,F7« 3«2A4« IHl » AZMP2840 

1B7 FORM AT ! 1H0.20X, SSHPROGRAM USES SUBROUTINE PART; CONSIDERS ONLY THAZMP2850 
lAT PORTION OF EACH VECTOR WITHIN THE CELL) A2MP2860 

199 FORMAT! 1H0,15X, 107HPROGRAM USES SUBROUTINE MIDI CONSIDERS WHOLE VAZMP2870 
lECTOR AS BEING WITHIN CELL IF ITS MIDPOINT FALLS IN THE CELL) AZMP2880 
190 FORMAT UH0,22HGR ID CELL NUMBERS ROW ,I4,9H, COLUMN ,I4,2H (,I4,5H AZMP2890 
KX< ,14, 2H; ,14, 5H <Y< ,I4,1H1) AZMP2900 

230 FORMAT !1H0,/,10H AZIMUTHS, 4X,16MCLASS LENGTH ! I N ,2 A4, 1 H ) ,37X , 1 9HAZMP2910 
INUMERICAL frequency, /, 2X, 8! 1H- ), 4X , 25! I H-) , 37X, 19 Cl M-) | AZMP2920 

205 FORHATCIH ,F5# I , 1 H- , F5# 1 , 2X,F8«2,57X ,14) AZMP2930 

216 FORMAT! 1H^.24X,IH>I AZMP2940 

223 FORMAT ! 1H+,87X, IH>) AZMP2950 

240 FORMAT! IH ,14X,7H-- ™ — ,58X, 4H -1 AZMP2960 

250 FORM AT < 1H0,5X, 6HTOTALS,2X,F8o 2 ,56X, I 5) AZMP2970 

265 FORM AT ! 1H*,28X ,9HEACH X = , F7« 2 , I X , 2 A4 ) AZMP2980 

295 FORMAT ( IH>,95X,9HEACM * = ,I4,6H UNITS) AZMP2990 

AZMP3000 
AZMP3010 


286 FORMAT! 1HF,24X,IH>,60A1) 

287 FORMAT! 1 H!> , 87X , I H>, 4 3A1 ) 


290 FORMAT (IH ,59H0NE OR MORE DENSITY HISTOGRAM CLASSES EXCEED LIMITS AZMP3020 


1 ALLOWED) 


AZMP3030 


295 FORMAT (IH ,61H0NE OR MORE FREQUENCY HISTOGRAM CLASSES EXCEED L I MIT A2MP30 40 
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c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 


WITHIN A GPIO 
IN THE CELL 


IS ALLOWED » 

297 FORM AT ( IHO, 13X ♦42HPANDDMLT DlSTRIflUTED OENSITT 0 
I f l2X«40HRAN0aMLT OISTRIBJTEO FREQ« DATA PROB,> = 
315 FORM AT C1M1*30X,71H*1TVPE* FOR PUNCHED OUTPUT NOT 
1E#A*44*J0B AB0RTE0*A444> 

330 WRITE! 1PR1NT.31SI 
400 STOP 
END 

SUBROUTINE PART 

SUBROUTINE *PAPT* DETERMINES IF A VECTOR FALLS 
CELL* AND IF SO* THE LENGTH OF THE PORTION WITH 

COMMON XI C200 0 >« YlC 20001 « X2< 2000 I « V2 ( 20 0 0 I • VECA ? 
10 ) *A (2000 )*0t 2000 If XM!0( 2 0001 * YM|0<2 000) .CLAMINI 
2AZLENC90KNAZFR0(90 I . XMI N . YMIN • XMA X t YMAX . FL AG * VL 
INTEGER VMINtXMlN»YHAX,XMAX*FLAG 
REAL MINLEN 

DATA I RE AO/5/* IPRINT/6/, M I NLEN /I .0/ 

DETERMINES IF WHOLE VECTDR IS WlTHlN GRID CELL 

ANY VECTOR < MINLEN (IN MM« > WILL NOT BE COUNTED 

IF(X1( Tl«GE«XMlN«ANOaXI<n*LE* XMAX« AND* Y 1 ( I IcGE* 
lY MAX* AND* X2C I ) »GE« XM IN* AND* X2( I >*LE* XMAX* AND* Y2( 
2t>*LE*YHAX) GO TO t 
GO TO 2 

1 FLAG =“l 
RETURN 

2 IFIAU n 3,20*3 

3 I F(A ( I l*LE*-573* *AND«X1( I ) «GE* XMAX) GO TO 422 
IF(X1( I)*EQ*XMAX*AND«A<n«GT*0* ) GO TO 20 
1F(X1( t )*6E«XMIN* AND«XH I ) «LT« XMAX, AND* Y 1 ( I I* GE» 

lYMAXl GO TO 100 

1 F(X2( X )• GE*XMXN«AND*X2( I ) *LT* XMAX* AND* Y2(I ) • GE» 
lYMAX) GO TO 200 
GO TO 300 

20 IFIXU I )*GE*XMTN*AND*XU t ) *LE* XM AX* AND* Y I ( I l*GE» 
lYMAX) GO TO 100 

IF(X2( t )*GE* XMIN*AN0*X2( X ) *LE* XMAX* AND* Y2 ( I )*GE* 
lYMAX) GO TO 200 
GO TO 300 


CALCU-^ATES VECLEN WHEN VECTOR ORIGIN IS WITHIN GRID CELL 

100 IFf A (I m 02* toi* 102 

101 VLENaX l( I )-XMlN 
FLAG*1 

RETURN 

102 IF (Xl< 1 )*EQ*XMlN*ANO*A(l )*GT*0* ) RETURN 
Y»A( I )4XMIN^B( I ) 

1F(Y*LE.VI(I )*ANO.Y. GE*Y2( I M GO TO 105 
GO TO 120 

105 I F(Y*GE*YMlN*AND*y*LE*VMAX ) GO TO 110 
GO TO 120 

110 ^/LEN=sSQRTC(XMIN-XlC I ) )*»2«^< V -YlCI))*42l 
GO TO 440 

120 Y sA( I l^XMAX^'Bl I ) 

!F(V*LE*Yl(t )*AND*Y*GE*Y2( IM GO TO 125 
GO TO 140 


AZMP3050 

AT A PPOB* = ,E-10^*4AZMP3 060 

«EL0*4) ' AZMP3070 

IN SPECIFIED RANGAZMP30B0 
AZMP3090 
AZMP3100 
AZMP3110 
AZMP3120 
PRTOOOlO 
PRT00020 
PPT00030 
PRT00040 
PRT00050 

M(2000) *VEClENC 200PRT00060 

90 ) *CLAMAX{90 ) , PRT00070 

EN,I PRTOOOeO 

PRT00090 
PRTOO 100 
PRTOOl 10 
PRT00120 
PRTOOl 30 
PRTOO I 40 
PR TOO 150 
PRTOO 160 

YM IN* AND* VI ( I )*LE* PRTOO 170 

! )*GE*YMI N* AND* Y2(PPT001B0 
PRT00190 
PRT00200 
PRT00210 
PRT00220 
PRT00230 
PRT00240 
PPT00250 

YMIN* AND* VI ( U*LT* PRTOO 260 
PRT00270 

YMIN* AND* Y2C 1 ) *LT*PRT002BO 
PRT 00290 
PRTOO 300 

YMIN* AND* Yl ( 1 ) • L T. PRTOO 3 10 
PRT00320 

YMIN* AND* Y2( I ) • LT* PRT 00 330 
PRTOO 340 
PRTOO 350 
PRT00360 
PPT00370 
PRTOO30O 
PRTOO 390 
PRT00400 
PRT00410 
PRT00420 
PPT00430 
PPT00440 
PRT00450 
PRT00460 
PRT00470 
PRT004B0 
PRT00490 
PRT00500 
PPT00510 
PRT00S20 
PRT00530 
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125 1F(Y •GF* YMIN • ANO. Y . LE. YMAX> GO TO ISO 
GO TO 140 

150 YLEIM=SQRT< <XMAX- XIC I n*|t24-( Y -Vlf! n**2> 

GO TO 440 

140 X-<YMIN-0( n 1/A{ 1 > 

IF(A(I I >141« 142» 142 

141 XFlXflGE^XlCl UAN0«X«LEaX2( in GO TO 145 
GO TO 170 

142 IF(X«LE«X1(1 )« AMO»X« GE«X2( I II GO TO 145 
GO TO 170 

145 IF (X «GE» XMIN * ANO * X »LE« XMAXI GO TO 165 
GO TO 170 

165 VLEN=SQRT(tX - X 1 ( I < YMI N- YU 1 1)44 2 I 

GO TO 440 

170 I F lA n I.GT^-STB* ) GO TO 175 
VLEN^Y It 1 l-YMIN 
GO TO 440 

175 IF< All )«LT«-'«0 ia« AND«A( I )*GT*0* > GO TO 180 
YLEN^XMAX«*Xt ( T ) 

GO TG 440 

130 !F<AC1 )*LT«0«.AND«A< I)«GT«*01d) GO TO 165 
VLENs=XlCI )-XMlN 
GO TO 440 

165 tfAlTEl TPRINT* 100 I XI C I ) « Y I C lU X2( I I • Y2C I I 
GO TO 422 

CALCULATES VECLEN mHEH VECTOR END IS WITHIN GRIO CELL 

200 IF<AU ) |202»201*202 

201 VLEN-XMAX-X2( J > 

GO TO 440 

202 Ya:A< I)4XMlN4eC I ) 

IFtV.LE.YK I J.ANO.Y. GE.Y2I I n GO TO 205 
GO TO 220 

205 1F(Y •GE* YWIN « AND* Y »LE* YHAX ) GO TO 210 
GO TO 220 

210 VLEN=SORTUXMIN-X2CI ) )4*24<Y -Y2(in442) 

GO TO 440 

220 Y-A( !»4XMAX+0{ I) 

IF(Y«LEoYl (I )• AND*Y« CE* Y2{ I n GO TO 225 
GO TO 240 

225 1FCY*6E*YNIN«AN0.V.LE*VMAX) GO TO 250 
GO TO 240 

250 VLEN=SQRT(<XMAX-X2< I n442U Y -Y2C1)>4*2) 

GO TO 440 

240 X = (YMA X^Bf I ) )/A< I ) 

1F1A< I )l24lt 242* 242 

241 I F(X«GE«X1< I )• AND«X«LE*X2C 1 ) I GO TO 245 
GO TO 270 

242 IF<X«LE«X1IX UANO«X*GE«X2( n > GO TO 245 
GO TO 270 

245 IF(X«GE<iXHIN«AND«X*LE*XMAX) GO TO 265 
GO TO 270 

265 VLeN=SQRTC<X - X2C 1 I 1*4 24< YMA X- Y2( I ) I *♦ 2 ) 

GO TO 440 

270 IF1A< 1 |*GT*-573* ) GO TO 275 
VLENi=YMAX-Y2< I I 
GO TO 440 

275 IF(A< I l«LT«^«O10«AND*A< I )«GToO« ) GO TO 280 
VLEN=X2( I I-XMIN 
GO TO 440 


PRT00540 
PRTOOS50 
PRT00560 
PRT00570 
PRT00580 
PRTQ0500 
PRT00600 
PRT00610 
PRT00620 
PRT00630 
PRT00640 
PR TOO 6 50 
PRT00660 
PPT006T0 
PRT006B0 
PRT00690 
PRT00700 
PRT00710 
PPT00T20 
PRT00730 
PRT00740 
PRT00750 
PRT00760 
PRT00770 
PRT00780 
PRT00790 
PRTOOeOO 
PRT 00810 
PRT00820 
PRT00830 
PRT00840 
PRT00650 
PRT00660 
PRT00870 
PRT00880 
PRT00890 
PRT00900 
PRT00910 
PRT00920 
PRT00930 
PRT00940 
PRT00950 
PRT00960 
PRT00970 
PRT00980 
PRT0Q990 
PRTOIOOO 
PRTOl 010 
PRTOl 020 
PRT01030 
PRT01040 
PRT01050 
PRT01060 
PRTOl 070 
PRT01080 
PRTOl 090 
PRTOl 100 
PRTOl 1 10 
PRTOl 120 
PRTOl 130 
PRTOl 140 
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c 

c 
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200 tF(A( I )«LT»0*«AND«A( ll.GT. .Oidl GO TO 205 
VLEN=XMAX«X2C I } 

GO TO 440 

205 WRITE(1PRINT*290) XS ( U • YM I I » X2C 1 I • Y2C t 1 
GO TO 422 

«PASS» DETERMINES WHETHER VECTOR PASSES THROUGH CELL 
AND IF SO« ITS LENGTH WITHIN THE CELL 

300 FLAGaO 
NCLC1=0 
NCLC 2=0 
NCLC 3=0 
NCLC 4=0 

IF(YH I )#LE«YMIN I RETURN 
IFIYU I I.LTaYMIN ) RETURN 
IFiY2C t U GE* YMAX ) RETURN 


PRTOl 150 
PRTOl 160 
PRTOl 170 
PRTOl 180 
PRTOl 190 
PRT01200 
PRT01210 
PRTOl 220 
PRT01230 
PRT01240 
PRTOl 250 
PRT01260 
PRT01270 
PRTO120O 
PRT01290 
PRTOl 300 
PRT01310 


IF<y 1C I UGE« VM INLAND* Yl( I ) «LE • YMAXik A NDa XI (I ) • LE» XM1N« AND* A( I ) * GT *OPRT01 320 
1*) RETURN PRT01330 

1F<Y1C 1 )*GE*YMIN*AND*Y1CI » *LE* YM AX* AND* X 1 ( 1 >*GE* XHAX* AND* A(I »*LT*0PRT01 340 
1* > RETURN PRT013S0 

IF|Y2( I I* GE. YMIN*AN0*Y2I I » *LE* YMAX* AND* X2 (I >*LE. XMIN* AND* A< 1 > * LT *OPRT 01 360 
1« > RETURN PPT01370 

IF<Y2< I I*GE*YMIN«AN0*Y2II ) *LE* YHAX« AND* X2(I )* GE» XM AX* AND* Alt U GT«OPRTOl 380 
1*1 RETURN PRTOl 390 

IFfYll I l*GE*YHAX*ANO*Xll 1 >*GE* XMAX*ANO* A( n*LT*0*l RETURN PRTOl 400 

IFIYII I »*GE* YMAX*AND«X1I I > *LE* XMl N* AND* A | II * GT« 0 . ) RETURN PRT01410 

IFCY 2( 1 )*LE* YMIN* AND* X2( T I *LE* XHTN* AND* A f I I * LT« 0* I RETURN PRTOl 420 

TF(Y2( I )*LE* YM1N*AND*X2( f I *GE* XMAX* AND* A ( II *GT*0*l RETURN PRT01430 

TF(Y2< t l*L£* YMIN* AND* X2( I I *LE* XM IN* AND* A ( I I * LT* 0** AND* A C I ) *GT* PRTOl 440 

1-573*1 RETURN PRTOl 450 

IFIACX l•^0*0*0*AND«Xlf II*GE*XMAX* AND*X2< I 9* LE*XMlN*ANDi Y1 <11 *GE* YMPRT01460 
1 t N*ANO«Y2( ! I*UT* YMAX 9 GO TO 301 PRT01470 

IF (A (I ),EQ*0*0«AN0*X2m»GE*XMAXI RETURN PRTOl 460 

IFCAI I )*EO*0*0*AND» XI ( I I*L£*XM!N I RETURN PRT01490 

IF(A<n*GT*-573* I GO TO 303 PRT01500 

1F(X1< 19* GE*XMIN*AND*XH ! I*LT« XMAX* AND* X2 ( I I « GE* XMIN* AN0*X2f I)*LT*PRTQ1 510 
1XMAX*AND* Yl< I l«GE*YMAX*AND*Y2( n*LE*YMlN*AND*An l»LE*-573* 9 GO TO PPT01S20 


1302 

PRTOl 530 


RETURN 

PRTOl 540 

301 

VLEN=XMAX-XMIN 

PRT01550 


GO TO 440 

PRTOl 560 

302 

VLEN=YMAX-YMIN 

PRTOl 570 


GO TO 440 

PRT01580 

303 

YXMI N= A< I 94XMlNTa( I 9 

PRTOl 590 


IF<YXM IN*LE* VI Cl >*ANO*YXMIN«GE*Y2< I n GO TO 305 

PPT01600 


GO TO 320 

PRTOl 610 

305 

TF<Y XM IN*GE* YMIN* AND*YXMIN*LE* YMAX 1 GO TO 310 

PRTOl 620 


GO TO 320 

PRT01630 

310 

NCLC 1=1 

PRTOl 640 

320 

YXMAX=A< ! 94XMAX>B<1 1 

PRTOl 650 


IF (YXMAX*LE* Y1 (1 l*AND*YXMAXcGE* Y2( I > 1 GO TO 325 

PRTOl 660 


GO TO 340 

PRTOl 670 

325 

IFIV XMAX* GE* YMIN*AND*YXMAX*LE*YMAX) GO TO 330 

PRTOl 680 


GO TO 340 

PRT01690 

330 

NCLC2=l 

PRTOl 700 


IF (NCLC1*GT*0«AND«NCLC2*GT*0I GO TO 335 

PRT01710 


GO TO 340 

PRT01720 

335 

VLEN=SORT< ( VXMAX-VXM IN 14112 + 1 XMI N-XMAX 2 1 

PRT01730 


GO TO 440 

PRT01740 

340 

XYMIN= < YM1N>»BC 1 1 I/AC 1 > 

PRTOl 750 
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c 

c 

c 

c 


IFCACI > )34l* 342*342 

341 1F< XVM tN« GE« XI (I )*AND*XYMIIN*LE*X2( t ) ) GO TO 345 
GO TO 300 

342 TF(KYMTN«LE* XUl UANO.XYmN«G£« X2C1 9 I GO TO 345 
345 tF(XYMtN<iGE«XMlM«AND«XYMINeLEft XHAX) GO TO 350 

GO TO 380 
350 NCLC3-1 

rFCNCLC2«GT« 0« AND«NCLC3* GT#09 GO TO 360 
GO TO 365 

360 VLE?^ = SQPT( < XM AX-XYMI N >♦* 24 ( YMl N-YXM AX I *4 2 ) 

GO TO 440 

365 tFfNCLCl«GT#0* AN0*MCLC3.GT«0I GO TO 370 
GO TO 380 

370 VLEN=SQRT< <XMIN-XVMIN9*42+< YXMIN- YMIN944 2 9 
GO TO 440 

3S0 XYMAXs ( YMAX- B( t 9 )/A< I 9 
IFCACI 9 9301* 362*362 

381 TFCXYMAX»GE.X1{I 9.AND*XYNAX«LE*X2CI>) GO TO 365 
GO TO 422 

382 IF(XYMAX«LE« XI ( I 9* ANO« XYMAX.GE* X2< t 9 9 GO TO 385 
60 TO 422 

365 IF(XYMAX* GE»XMIN«ANO*XYMAX*L£« X»4AX9 GO TO 300 
FU AG=0 
RETU RN 
390 NCLC4-1 

IF<^CLC3#GT*0« AN0«NCLC4*GT«0| GO TO 395 
GO TO 400 

395 VLEH=sSQRT( (XYMAX-XVM IN9 ^4 24( YMAX*- YMI N 9 44 2 9 
GO TO 440 

400 TF(9CLC2.GT«0#ANO«NCLC4.GT*0> GO TO 405 
GO TO 410 

405 VLEN=SQRT( (XMAX-XYMAK94424(YMAX-VXMAX)*429 
GO TO 440 

410 IF (NCLC1«GT«0« AND«NCUC4« GT«0 9 GO TO 426 

IFiNCLCloGE. l*DR*NCLC2.GE* !• OR# NCLC3# GEo NCL 

420 WP1TE< IPRIMT*430 ) NCLCl ♦ NCLC2 • NCLC3 * NCLC4 ,X 1 ( I 9 • 
422 FLAG-0 
RETURN 

425 VLENs:50PT< {XYMAX-XM1N9A4 24C YMAX-YXMIN94429 
440 FLAG=1 

IF<VLEN#LT«MINLEN) FLAG=0 
RETURN 

190 FORM AT (lHO*70H4<^44*1r4*E^ ^0R44444444* VECTOR OPIG 
ITOR DELIMITERS ARE* 4 € 5X» F7* 1 9 * // * IH .50H44PRPOR 
2FOLL OWING PRINTED CELU44 9 
290 FORMAT ( lHO*67H44*4*444ERROR4fcA^A»4*w VECTOR END 
1 DELIMITERS APE * 4 ( 5X * F7» I 9 * // * 1 H ,60H44ERROR meS 
2LOWING PRINTED CELL449 

430 FORM AT C lH0*63H4*4444EPP09A*k4A44 SUBROUTINE VECTO 
lNCLCl-4 ARE* 41 5 *22H VECTOR DELIMITERS ARE*4F7. 

END 

SUBROUTINE MIO 

SUBROUTINE 'MID* DETERMINES IF VECTOR MIDPOINT F 
CELL 

COMMON XI 12000}* Yl( 20009 * X2( 20 00 I » Y2< 200 O 9 * VE C AZ 
10 > *4(2000 9*B(200 09*XM!D< 2000 9 *YM!0( 2000) *CLAMlN( 
2AZLENC 90 ) «NAZFRO C 90 9 * XMIN* YMIN* XMAX ■ YMAX *FL AG*VL 
INTEGER XMlN * YMlN*XMAXvYMAX*FLAG 

TF(XMI0( I )*GE«XMlN«AND»XMIO(I 9 *L T* XMAX« AND* YMtOC 


PRT01760 
PPT01770 
PRT01760 
PRT01790 
PRTOl 000 
PPTOieiO 
PRTOt 020 
PR TO 18 30 
PRTOl 840 
PRTOl 050 
PRTOl 060 
PRTOl 070 
PRT018S0 
PRT01890 
PRT01900 
PRT01910 
PRT01920 
PRT01930 
PRT01940 
PRT01950 
PRTOl 960 
PRTOl 970 
PRTOl 980 
PRT01990 
PPT02000 
PRT02010 
PPT02020 
PRT02030 
PRT02040 
PRT02050 
PRT02060 
PRT02070 
PPT02080 
PRT02090 

C4*GE*19 GO TO 420PPT02100 

Yl (1 9 *X2C I 9* Y2C1 9 PRT02110 
PRT02120 
PRT021 30 
PRT02140 
PRT02150 
PRT02160 
PPT02170 

IN SUBROUTINE* VECPRT021B0 

MESSAGE REFERS TO PRT02190 
PRT02200 

SUBROUT t NE* VECTORPRT022 10 

SAGE REFERS TO F0LPRT02220 
PRT02230 

R-PASS; VALUES OF PRT02240 

19 PRT02250 

PRT02260 
MIDOOOlO 
MT 000020 

ALLS WITHIN GRID MI000030 
MI D00040 
MIOOOOSO 

M(20009 *VEC_EN( 200MID00060 

90I*CLAMAX(90 9 * MI 000070 
EN*i Miooooeo 

M1D00090 

1 )«GE»YM1 N* AND* YMIMfDOOl 00 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 


10 < n •LT.yWAX > <^3 TO 1 
FLAG^O 
RETJ RN 
I FLAG-*1 
RETURN 
END 

FUNCTION PRBCHI fCHlSO«IOF> 

^IRITTEN BY H*D* KNO 0 LE & F. YATES aOPOEN# THE PENNSYLVANIA STATE 
UNIVERSITY* 1966 

THIS FUNCTION COMPUTES BY THE APPROXIMATIONS ON PAGE 9A1 OF 
^HANDBOOK OF MATHEMATICAL FUNCT IONS** *U« S* OEPT* OF COMMER CE* 1964 « 
GIVEN A VALUE OF CHI-SQUARE AND ITS DEGREES OF -REEDOM* FUNCTION 
PRBCHI COMPUTES THE PROBABILITY OF A GREATER VALUE OF CHI-SQUARE* 
THE Z( ARGUMENT) FUNCTION IS COMPUTED BY FORMULA 26*2«1* P# 931# 

INTEGER TEST 

ALL PE ALA a ARGUMENTS CHANGED TO DOUBLE PRECISION BY M* POO WVSOCK I • 
DOUBLE PRECISION OSQRT*DE XP* AR G. SCHI SQ* XPLE VL 
double precision Q*R »S#T*U. V* V9*PPOa,S2Pl *2006*APPR0X 
DATA S2PI/2*5066262000/ 

QC ARC) sf OEXP C-ARGAARGAO* S| /a* 506620 2D00 ) * C T A C Oa 3 1 9361 5000 ♦TA 
U-0» 35 6 5636000 +T AC 1« 761 476000+TA 1. 821 2563 00«M 1 * 330274000AT )) » ) | 

X PI,- 2e 576 235 9600 0 
PRBCHls0«0 

XF<CHISQ«LT«0«0) RETURN 
tFUDF«LE*0) RETURN 
100 SCHISQ^CHISQ 
S=lc 0 
V== IDF 

V 9*2 • 0 /FL O AT C 9 A I DF ) 

U*-SCHISQAOft 3 

SCHI S0=DS0RTCSCHISQ) 

IF (0ABSCU)*LT«174*6t GO TO llO 

174*6 IS THE largest ARGUMENT THAT EXP WILL TAKE# 

PROSsO#0 
GO TO 240 

CHECK FOR DEGREES OF FREEDOM GREATER THAN 100 OR GREATER THAN 30 

no IF CIDF#GT«100» GO TO 200 
IF (1DF#GT#30I GO TO 170 

DEGREES OF FREEDOM LESS THAN OR EQUAL TO 30 

PR OB *0*0 
TEST=MODt lOF* 2 > 

IF (TEST#NE#0) GO TO 140 

EVEN DEGREES OF FREEDOM AA LESS THAN OR EQUAL TO 30 AA FORMULA 
26#4«5* PAGE 941 

1RANG£=( IDF-2)/^2 

IF ( IRANGE.EQ# 0) GO TO 130 

DO 120 I*1*IPANGE 

I«=I +I 

S*SA TP 

120 PR08*PR0B+SCMISQAAIR/S 


MiDOOl 10 
MID00120 
Mi DOO 1 30 
MIDOO 1 40 
MI000150 
MX DOO 1 60 
CHIOOOX 0 
CHI00020 
CH100030 
CH100040 
CHX00050 
CHI00060 
CHI00070 
CHI 00080 
CHI00090 
CHIOO 100 
CHIOOllO 
CHI 00 120 
CHI00130 
CHI00140 
CHtOOlSO 
CHI00160 
CHIOO 170 
ICHIOOieo 
CHIOO 190 
CHI00200 
CH100210 
CHI00220 
CMl 00230 
CHI 00240 
CHI00250 
CH1Q0260 
CMI00270 
CH100280 
CHI 00290 
CHIOO 300 
CHI00310 
CHX00320 
CHI00330 
CHIOO 340 
CHX00350 
CHIOO 360 
CHI00370 
CHI 00380 
CHI 00390 
CHI 00400 
CH100410 
CMI00420 
CHI00430 
CHI00440 
CHI00450 
CHI 00460 
CH1004TO 
CHI00460 
CHI00490 
CHIOO 500 
CH100510 
CHI00520 
CH 100530 
CHI00540 
CH100550 
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130 


PR08-DEXP<U)*( l*0+PROB) 
GO TO 230 


C 

c 

c 

c 


000 DEGREES OP FREEDOM 
PAGE 041 


♦♦ CESS THAN OP EQUAL TO 20 ♦♦ FORMULA 


140 


150 

160 


IRANGE^C lOF-l 1/2 
IF ( IPANGE.£Q«0) GO TO 160 
DO 150 I = l» LBANGE 

IR-I ♦!-! 

S=S> IR 

PROBSPROB4SCH1S044IR/S 

T=l. 0/ ( 1 .040 *2 31 64100004 SC HI SOI 

PROQs2«04 CQCSCHISQ) ) 42.041 DEXPtU)/S2PI 14PR0B 

GO TO 230 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


****«»4»*4 GREATER THAN 30 DEGREES OF FREEDOM *4444*444* 

AH APPROXIMATE VALUE OF CHISO IS FIRST COMPUTED THEN COMPARED WITH 
THE GIVEN CHISQe IF THE APOrqx* VALUE IS GREATER THAN THE GIVEN 
VALUE, 0<CHISQ,IDFI IS RETURNED AS •995. 


CH100560 
CHT00570 
CHI 00580 
CHI 005 90 
CHI00600 
CHI00610 
CH100620 
CHI 00630 
CHI00640 
CH100650 
CHI 00660 
CHI00670 
CHIOO60O 
CH100690 
CHI00700 
CHI00710 
CHI00720 
CHf00730 
CHI00740 
CHI00750 


A44444 44*44«<***44ir *444 44444 4 44 4 444*4444*4 44^4 4*4*4 4444444*4444 44 44 **CHl 00760 
FOR GREATER THAN 30 AND LESS THAN OR EQUAL TO 100 DEGREES OF FREEDCMCS 1 00770 
THE APPROX. VALUE OF CHXS3 AT THE .995 LEVEL IS COMPUTED BY FORMULA CHI00780 
26.4.17, PAGE 941. THE SIGN OF X«PI IN THE FQPMULA WAS CHANGED 
FROM 4 TO - TO ALLOW COMPUTATION OF CHlSQ AT THE .995 LEVEL RATHEP 
THAN THE .005 LEVEL AS IS THE CASE WHEN THE SIGN IS 


170 


160 

190 


APRO Xs C (1.0‘-V9-XPL4OSQRTC V9n4431*V 
IF { APROX.LE.CHI SQl GO TO 180 
GO TO 210 

V=<( CHISQ/VI440. 33333333000-< 1.0-V9I 1/D$QRT(V9) 
T = l. 0/(1.040.2 3164190004V> 

PROBaCHVl 
GO TO 230 


GREATER THAN 100 DEGREES 0= FREEDOM. 
IS COMPUTED BY FORMULA 26.4.16, PAGE 
CHANCED FOR THE SAME REASDM AS ABOVE 


THE APPROX. 
941, THE SIGN 


VALUE OF CHISO 
OF X(PI WAS 


200 APR0X=U-XPL4DS0RT(V4V-1 .0 >)4*2>40.5 

IF ( APROX.LE.CHtSQI GO T3 220 
210 PRO0S4O.995 
GO TO 240 

220 VsOSORTI2.00 04CHISQ)'-OSQRT<2.O4V- 1.0 ) 
GO TO 190 

230 IF ( PROB.GT. 0*995 > GO TO 210 
240 PPBCHlsPROB 
RETURN 
END 


CHI00T90 
CHI00800 
CHlOOeiO 
CHI00820 
CHIOO03O 
CHI00840 
CHI 00650 
CHI00860 
CMI00870 

CHiooeeo 

CH100690 
CHI0090Q 
CHI00910 
CHI00920 
CHI00930 
CHI00940 
CHI009S0 
CHI00960 
CHI00970 
CH100980 
CHIQ0990 
CHIOIOOO 
CHlOlOlO 
CHI01020 
CHlOl 030 
CHIO1Q40 
CHIOIOSO 
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ROSE DIAGRAM PROGRAM 


PRBCS)ING^ 


page blank not FIL"KD 



c 

c 

c 

c 

c 

r. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


♦ POSEO 0 10 

RDSE DIAGRAM PROGRAM RD55E0020 

PH5;E0 030 
ROSE0040 

THE PROGRAM WAS WRITTEN! BY MELVIN H« PODWYSOCKI OF THE GEO- ROSE0050 

SCIENCES DEPT*# THE PENNSYLVANIA STATE UNIVERSITY. MAY. 1972 FORRDSE0060 
USE ON AN I8M 360/67 COMPUTER ALONG WITH THE PENN STATE COMPU- ROSE0070 
TAT ION CENTER «QDGS** GRAPHICS PACKAGE. THE PROGRAM WAS MODIF I EOROSEO 0 80 
IN MAY. 1973 TO USE STANDARD CALCOMP SOFTWARE SUBROUTINES WITH ROSEO09O 
THE CALCOMP 700 DRUM PLOTTER PACKAGE ON THE IBM 360 SERIES POSEO 100 

COMPUTERS HAVING THE EQUIVALENT OF 130K BYTES STORAGE* ROSEOllO 

ROSEOl 20 

PROQPAM PRODUCES ROSE DIAGRAMS SUITABLE FOR MAP OVERLAY FOR ROSEOl 30 
ANY NUMBER OF MAP GRID CELLS. PLOTTING THE ROSE A'*' THE MIDPOINT ROSE0140 
□ p XHE cell along WITH THE SUM TOTAL OF THE VALUES COMPRISING RDSEOISO 

THE ROSE* SCALING FACT3RS ARE INClOdED WITHIN THE PROGRAM, PQSE0160 

DATA INPUT IS GENERATED BY THE ••AZMAP** PROGRAM AND CONSISTS OF X R0SE0170 

AND Y MIDPOINTS (XMIO S YMIDI OF E ACH GRID CELL AND THE COMPO- RDSEOiaO 
NENTS OF THE ROSE DIAGRAM (AZLEN!o UP TO 90 AZIMUTH CLASSES MAY RDSE0190 
BE USED BETWEEN 270 THRU 0 TO 90 DEGREES* COORDINATES MUST BE RDSE0200 

READ IN MMo . X 1 S ♦ TO THE RIGHT AND Y IS + DOWNWARD* NORTH IS RnSE0210 

ASSUMED PARALLEL TO THE V AXIS* CONTROL CARDS ARE READ FROM THE ROSE022O 
CARD READER WHILE DATA CAROS APE READ FROM ANY UNIT DECLARED BY R05E0230 
•ITAPEl* IN STATEMENT RDS 660* OUTPUT IS GENERATED AS PER THE R0SE0240 

CALCOMP PACKAGE AT EACH INDIVIDUAL INSTALLATION AND CONSISTS OF R0SE0250 
A PLOT OF X S Y AXES AND ROSE DIAGRAMS* THE TOTAL VALUE OF UNI TSROSEO 2 60 
COMPRISING EACH ROSE DIAGRAM IS ALSO PLOTTED* R0SE0270 

ROSE028O 

ALL NUMERIC INPUT DATA IS RIGHT JUSTIFIED! «I« INDICATES I NTEGER. POSE0290 
tipt. indicates FLOATING >DINT AND "A** INDICATES CHARACTER FORMAT; ROSE0300 
PPECEEDING NUMBERS INDICATES COLUMNS USED = OR EACH PARAMETER* ROSEO 3 1 0 


ROSE0320 

♦ CONTROL CARD 1 — — PARAMETER CARD R0SE0330 

XMN*M1nIMUM X VALUE FOR »LOT IN MM* <F7, 2,4 1-7) RDSE0340 

YMNaMlNIMUM Y VALUE FOP ^LOT IN MM* < F7. 2 , 4 B- 1 41 RDSE0350 

XMX=MAXIMUM X VALUE FOR PLOT IN MM* CF7* 2 . 4 1 5-21 ) ROSE0360 

YMXxMAXIMUM Y VALUE FOR PLOT IN MM* < F7* 2 , 4 22-28 ) PD5E0370 

note: THE ABOVE 4 VALUES GOVERN THE X AND Y AXES LABELING* ROSEO 380 

WHICH WILL 0E INCREMENTED BY THE VALUE • SC • C ROS 6501* SO ROSEO 390 

THAT WHEN ’FACT* = 1* 1 INCH = *SC* MM* ROSE0400 

F ACT=MULTIPL ICATION FACTOR FOR ALL CARTESIAN COORDINATE DATA; AL- ROSE0410 
LOWS SCALING LARGE DIMENSIONED MAPS TO SIZE ACCOMMODATED BY R0SE0420 
PLOTTER* MAY BE <* = OR > 1 • BUT NOT < 0* < F 7* 2. 429-35 I R0SF0430 

RSSZ-MULT tPL ICAT ION FACTOR FOR *A2LEN» OF ROSE DIAGRAM* SCALES POSE0440 
ROSE DIAGRAMS SO THAT THEY DO NOT OVERLAP EACH OTHER OR EXCEED ROSE04S0 
BOUNDS OF PLOT AND IS EMPERICALLY DEVELOPED* C F7, 2 . *36-42 ) R0SE0460 

note: ALL *AZLEN* VA_UES ARE SCALED BY 'PACT* AS WELL AS •SC* ROSE0470 
NCELLttNUMBER OF MAP GRID CELLS <14,443-46) R0SE0480 

NCLASS=NUMBER OF AZIMUTH CLASSES <1*E« 180 OEGPEES/10 DEGREE CLASSROSE0490 

INCREMENT = 18 CLASSES) (12.447-48) ROSEO 500 

PGS-Y AXIS PAGESIZE OF CALCOMP PLOTTER <IN INCHES. t*E* 12**30*. R0SE0510 

ETC*) CF5* 0. 449- S3) RDSE0520 

«4t 4*^ CONTROL CARD 2— —--TITLE CARD R0SE0530 

TITLE WILL BE PLACED AT TOP OF PLOTTED OUTPUT* IF CENTERING IS R0SE0540 

DESIRED. IT MUST BE pJNCHED SYMMETRICALLY ABOUT *40 OF THE POSE0550 

TITLE CARD (20A4. 41-80) : ROSE0560 

**f CONTROL CARO 3 ^-INPUT FORMAT CARD ROSEOS70 


FORMAT MUST 0F ENCLOSED IN PAPANTHESES AND BEGIN IN 41 * SEQUENCE ROSE0580 
MUST be: XMlD. VMIO. AND THE NUMBER OF • AZLEN • CLASSES AS SPE- P0SE0590 


ClFieO BY 'NCLASS* (I*E* < 2F6* 2. 8F 6* 2 / 10 F 6. 2) FOP 'NCLASS* = R0SE0600 


ritiiSOiiDiiSlXi PAGE BLANK NOT PILMELl 
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c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


l«l (20A4.1 0O> ROSF0610 

****„*<r*OATA CARDS-- RDSF0620 

READ FROM ANY UNIT DECLARED BY 'ITAPEt** MUST BE IN FORMAT SPEC 1-ROS E0630 
FIED BY CONTROL C^RO 3. R0SE06A0 

ROSF0650 
ROSE0660 
POSE0670 
ROSEO60O 
PPSF0690 
ROSeOTOO 
ROSEOTIO 
R0SE0720 
POSE0730 
ROSE0740 
R0SE0750 
POSF0760 
RnSEOTTO 
ROSE0780 
R0SE07R0 
ROSE0800 
ROSEOBIO 
ROSE0820 
ROSEO03O 
P0SE0840 
RnsEoeso 
P0SE0660 
RDSE0670 
ROSEOftSO 
ROSF0890 

CALL AXIS(0. ♦O. W3HX AXIS OF MAP • I 3 . C XMX- XMN • /SC • 0. , XMN , 50. * 1 Oo > ROSE0900 

_ ROSE0910 

PLOT NORTH ARROW PARALLEL TO Y AXIS (SIDE OF PAGE) RQSE0920 

ROSE0930 
R0SE0940 
ROSE0950 
PnSE0960 
ROSEO 970 
«nSE0980 
POSF0990 
RGSEIOOO 
POSElOlO 
ROSEl 020 
R0SE1030 
POSEl 040 
ROSEIOSO 
pnsEioeo 

^ RnSElOTO 

CALL AXISiO# tO,, 13HY AXIS OF M AP 1 3 * ( YMX-Y MN » /S C » -90« . YMN, 5 0* • I Oo ROSEl 080 

POSEl 090 


D1 WENS ION AZLENt 1 80) . AZMI D( I 80 ) « AZ X ( 1 80 ) • AZY( 180 ) t ABXU 80 ) » 
1 A8Y( ISO )• A2R A0< 1 80 ) 

dimension buffer (SO 00) .FMTROt 20 I^TITLE < 20) 

DATA PRAD/U745329E- 2/f I R EAO/5/* SC/50* /« HT25/* 25/, HTI /, 1/, 

1 lTAPEl/6/ 

CALL PLOTSIBUFFERtSOOO I 
READ CONTROL CARDS 

READ C I READ, 5 I XMN,YMN,XMX, YM X, FACT , R SS2 , NCE LL ,NC LASS . PGS 
READ<inEAD*10) ( T ITL E<L ) • L -I , 2 0) 

RE AD (I PFA0,1 0 ) ( FMTPD(L) ,L = t ,20 I 

MOVE PEN TO TOP OF PAGE (ASSUMES PEN SET TO PAGE BOTTOM) 

PGSZxPGS-Oe I 54PGS 
CALL PLOT (0* ,PGS2,23 ) 

FACTOR THE FOLLOWING PLOTTING SUBROUTINES 
CALL FACTOR(FACT) \ 

PLOT X AXIS PARALLEL TO TOP OF PAGE 


ANRT X=: xmx/sc*i«o 
ANR r=XWX/SC^I*S 
call PL0T(ANRTX,-3* • 31 
CALL PLOT(ANRTX,-l« 5,2) 

CALL PLOT ( AMRT 2* • 2 ) 

CALL SYMBOL ( ANRT X-HT 25/4, , - 1 * 25, HT25* IHN , 0» , 1 ) 
PLOT TITLE ABOVE X AXIS 


TSTRT^C <XMX-XMN)/{SC42o))-0* 5«'80*n:HT25 
CALL SYMBOL! TSTRT, I * 25,HT 2 5, T I TLE , Oe ,80) 

PLOT Y AXIS PARALLEL TO SIDE OF PAGE 


1 ) 

NC=27NCL ASS 


READ DATA FROM UNIT ITAPEl , ONE ROSE DIAGRAM AT A TI 
DO 200 K=1,NCELL 

READ (ITAPEl, FMTRO ) XMTD, YH 10, ( A ZLENI M ) , M = I , NCLA5S) 
T0TLN=0. 

PLOT EACH ROSE DIAGRAM 

DO 140 I=fl,NC 
NCl=NCLASS^'I 


ME 


ROSEl too 
Rosei 110 

POSEl 120 
ROSEl 130 
ROSEl 140 
ROSEl 150 
POSEl 160 
ROSE! 170 
ROSFl 180 
ROSEl 190 
POSEl 200 
R0SE12I0 
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onn n n n nno onr» onn 


15 

^0 


60 

flO 


100 


1 20 
140 


20 0 


5 

10 


IF(I-WCLASS» IS, 15.100 
t F(T-1 I 40.40,60 
ySCALF=“YMlD/SC 
XMID=XMID/SC 

PLOT CENTEPPOINT OF EACH ROSE DIAGRAM 

CALL symbol ( XMID“HT l/4«, YSCALE-HTl/S^.HTaS^RSSZ, 1H4.0* • 1 } 

IF( I-NCLASS > 00,80,100 

AZMI0( I )a270.- 130*/NC4lt«180. /NCLASS 

I F(AZMT0< I UGT.360* ) AZMIOU ) = A ZMIO ( T )«360, 

A ZMl D( NCI )=AZM 1D( r 14 180:, 

lF(AZMIOCNCl UGT.360* > AZMIDCNCI )*AZMlD<NCl >-360« 

AZL£N( NCI )5fA2LEN< T ) 

T0TLN*T0TLN4 AZLENU > / 

TF(AZM 10( I )* GF*0**AND, A2M!D< I ULToOGo )' AZMI 01 I I = 90 ® - AZM 1 0 C I » 

TF ( AZMlDf I )» GF«O0««AN0«AZM ID Cl ULT 9 366« ) AZMIDC I 1 -450 • - AZ MID C I ) 
AZRADC I ) = AZM IOC I M«PRA0 

SCALE AZLEN BY RSSZ AND SC 

ABXC II = C AZLENC I )*COS C AZRADC I I >4PSSZ)/SC 
AZXC I IsABXC I )4XM ID 

ABYC U = ( AZLENC n *SIN C AZRADC 1 ) )4R$SZ|/SC 
AZYC n=:ABYCI |4YSCALE 

PLOT EACH AZLEN 

KEYS 2 

IFCI«EQ«1 I KEY=3 

CALL PLOTCAZXtlJ ,AZY(II,<EY) 

IFCI-NCl 140,120,120 
CALL PLOT CAZXC II ,AZV< 1 1.2 I 
CONT INUE 

S TRT =A ZXC NCL ASS I 40o 5ARSSZ 

PLOT AZLEN SUM FOR EACH ROSE DIAGRAM 

CALL NUM0ERCSTRT, AZYCNCLASS ). 2« ♦HT 25*RSS Z ,T OTLN, 0« ,1 1 
CONTINUE 

TERMINATE PLOT 

CALL PLOT (XMX/SC+5# , C VMM - V MX > /SC .999 ) 

FORMAT C6FT* 2. 14. I2,F5«0> 

FORM AT { 20 A4) 

STOP 

END 


ROSEl 220 
ROSE! 230 
ROSEl 240 
ROSEl 250 
ROSEl 260 
ROSEl 270 
ROSEl 260 
ROSEl 290 
ROSEl 300 
ROSE1310 
ROSEl 320 
ROSE I 3 30 
ROSEl 340 
ROSEl 350 
ROSEl 360 
ROSEl 370 
ROSEl 360 
ROSEl 390 
ROSE1400 
ROSE1410 
ROSEl 420 
ROSEl 430 
ROSEl 440 
ROSEl 450 
RDSE1460 
ROSEl 470 
ROSEl 460 
ROSEl 490 
ROSEl 500 
ROSEl 510 
ROSEl 520 
ROSEl 530 
ROSEl 540 
ROSEl 550 
ROSEl 560 
ROSEl 570 
POSE! 560 
ROSEl 590 
ROSE1600 
ROSE1610 
ROSEl 620 
P0SE163O 
ROSEl 640 
R0SE1650 
R0SE1660 
ROSEl 670 
P0SE1680 
ROSEl 690 


NASA-GSFC 
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